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ABSTRACT 

The observed two-point angular correlation function, of mJy radio sources 

exhibits the puzzling feature of a power-law behaviour up to very large (~ 10°) an- 
gular scales which cannot be accounted for in the standard hierarchical clustering 
scenario for any realistic redshift distribution of such sources. After having discarded 
the possibility that the signal can be explained by a high density local {z < 0.1) 
source population, we find no alternatives to assuming that - at variance with all 
the other extragalactic populations studied so far, and in particular with optically 
selected quasars - radio sources responsible for the large-scale clustering signal were 
increasingly less clustered with increasing look-back time, up to at least z ~ 1. The 
data are accurately accounted for in terms of a bias function which decreases with 
increasing redshift, mirroring the evolution with cosmic time of the characteristic halo 
mass, M^,, entering the non linear regime. In the framework of the 'concordance cos- 
mology', the effective halo mass controlling the bias parameter is found to decrease 
from about 10^^ Mq/H at z ~ to the value appropriate for optically selected quasars, 
~ 10^^ Mq/H, at z ~ 1.5. This suggests that, in the redshift range probed by the data, 
the clustering evolution of radio sources is ruled by the growth of large-scale structure, 
and that they are associated with the densest environments virializing at any cosmic 
epoch. The data provide only loose constraints on radio source clustering at z > 1 
so we cannot rule out the possibility that at these redshifts the clustering evolution 
of radio sources enters a different regime, perhaps similar to that found for optically 
selected quasars. The dependence of the large-scale shape of ■w{9) on cosmological 
parameters is also discussed. 
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1 INTRODUCTION 

Extragalactic radio sources are well suited to probe the large 
scale structure of the Universe since they are detected over 
large cosmological distances (up to 2 ~ 6), are unaffected by 
dust extinction, and can thus provide an unbiased sampling 
of volumes larger than those usually probed by optical sur- 
veys. On the other hand, their 3D-space distribution can be 
recovered only in the very local Universe {z £ 0.1; see Pea- 
cock & Nicholson 1991; Magliocchetti et al. 2004) because 
the majority of radio galaxies detected in the available large 
area surveys, carried out at low frequencies, have very faint 
optical counterparts, so that measurements of their redshifts 
are a difficult task. As a result, only the angular clustering 
can be measured for the radio source population. Interest- 
ingly, high-frequency surveys have much higher identifica- 



tion rates (Ricci et al. 2004), suggesting that this difffculty 
may be overcome when such surveys will cover sufffciently 
large areas. 

Even the detection of clustering in the 2D distribution 
of radio sources proved to be extremely difficult (see Webster 
1976, Seldner & Peebles 1981, Shaver & Pierre 1989) since, 
when projected onto the sky, the space correlation is sig- 
nificantly diluted because of the broad redshift distribution 
of radio sources. It was only with the advent of deep radio 
surveys covering large areas of the sky, such as the FIRST 
(Becker, White & Helfand 1995), WENSS (Rengelink et al. 
1997), NVSS (Condon et al. 1998), and SUMSS (Bock, Large 
& Sadler 1999), that the angular clustering of this class of 
objects has been detected with a high statistical significance 
down to flux density limits of few mJy (see Cress et al. 1997 
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and Magliocchetti et al. 1998, 1999 for the FIRST survey; 
Blake & Wall 2002a,b and Overzier et al. 2003 for NVSS; 
Rengelink & Rottgering 1999 for WENSS and Blake, Mauch 
& Sadler 2004 for SUMSS). Amongst all the above surveys, 
NVSS is characterized by the most extensive sky coverage 
and can thus take advantage of high statistics despite its 
relatively high flux Umit (~ 3mJy vs ~ ImJy of FIRST). 
The two-point angular correlation function, w{0), measured 
for NVSS sources brighter than 10 mjy is well described by 
a power- law, extending from 0.1 degrees up to scales of 
almost 10 degrees. A signal of comparable amplitude and 
shape was detected in the FIRST survey at the same flux 
density limit, on scales of up to 2-3 degrees (see e.g. Maglioc- 
chetti et al. 1999), while at larger angular separations any 
positive clustering signal - if present - was hidden by the 
noise. 

Most analyses of radio source clustering performed so 
far (see e.g. Blake & Wall 2002a,b; Overzier et al. 2003) as- 
sume a two-point spatial correlation function of the form 
^rg{r) = (r/ro) The power- law shape is in fact preserved 
when projected onto the sky (see Limber 1953), so that the 
observed behaviour of the angular correlation is well recov- 
ered. The studies summarized by Overzier et al. (2003) typi- 
cally found correlation lengths ro in the range 5-15 h^^ Mpc. 
This large range may reflect on one hand real differences 
in the correlation properties of radio sources of different 
classes/luminosities, and, on the other hand, the large un- 
certainties on both the redshift distribution of the sources 
and the time-evolution of clustering, which are necessary 
ingredients to estimate $rg(r) from the observed w{0). 

Overzier et al. (2003) found that the observed cluster- 
ing of powerful radio sources is consistent either with an 
essentially redshift independent cornoving correlation length 
ro = 14 ± 3h~^ Mpc, close to that measured for extremely 
red objects (EROs) at 2: ~ 1 (Daddi et al. 2001) or with 
a galaxy conservation model, whereby the clustering evo- 
lution follows the cosmological growth of density perturba- 
tions. In the latter case, the present-day value of ro for the 
most powerful radio sources is comparable to that of local 
rich clusters. 

A deeper examination of the power-law behaviour of the 
angular two-point correlation function up to scales of the 
order of ~ 10° highlights some interesting issues. In fact, 
within the Cold Dark Matter paradigm of structure forma- 
tion, the spatial correlation function of matter displays a 
sharp cut-off around a comoving radius of r ~ 100 Mpc* 
(see e.g. Matsubara 2004, Fig. 1), which, at the average red- 
shift for radio sources < z >~ 1, corresponds to angular 
separations of only a few (~ 1° — 2°) degrees, in clear con- 
trast with the observations. This opens the question of how 
to reconcile the clustering properties of these sources with 
the standard scenario of structure formation. Magliocchetti 
et al. (1999) claim that the large-scale positive tail of the 
angular correlation function w(0) can be reproduced by al- 
lowing for a suitable choice of the time-evolution of the bias 
parameter, characterizing the way radio galaxies trace the 
underlying mass distribution. Although promising, this ap- 



* We assume a flat universe with a cosmological constant and 
Ho.m = 0.27, Clb = 0.045, as = 0.9, n = 1, and /i = 0.72, in 
agreement with the first- year WMAP results (Spergel et al. 2003). 



proach suffers from a number of limitations due to both 
theoretical modelling and the quality of data then available. 

The aim of the present work is therefore to investigate 
in better detail and to provide a self-consistent explanation 
of the puzzling behaviour of the angular correlation function 
of radio sources, especially on large angular scales. We will 
concentrate on the results from the NVSS survey (Blake & 
Wall 2002b; Overzier et al. 2003) as, thanks to the extremely 
good statistics, a clear detection of a positive clustering sig- 
nal was obtained up to ~ 10°. We will exploit the available 
spectroscopic information on local radio sources to limit the 
uncertainties on their redshift distribution, and will mainly 
focus on the time evolution of their clustering properties via 
the bias parameter. In this way we will derive interesting 
constraints on the typical mass of dark matter halos hosting 
the population of radio sources. We will also investigate the 
dependence of the predicted angular correlation function on 
cosmological parameters. 

The layout of the paper is as follows. A short descrip- 
tion of the NVSS survey is presented in Section 2. Section 3 
illustrates the adopted model for the redshift distribution of 
mJy radio sources, while in Section 4 we provide the formal- 
ism for the two-point angular correlation function. Results 
and discussions are presented in Section 5. In Section 6 we 
summarize our main conclusions. 



2 THE NVSS SURVEY 

The NVSS (NRAO VLA Sky Survey, Condon et al. 1998) is 

the largest radio survey that currently exists at 1.4 GHz. It 
was constructed between 1993 and 1998 and covers ~ 10.3 sr 
of the sky north of 5 = —40°. The survey was performed 
with the VLA in the D configuration and the full width 
at half-maximum (FWHM) of the synthesized beam is 45 
arcscc. The source catalogue contains 1.8x10*^ sources and 
it is claimed to be 99 per cent complete above the integrated 
flux density Si.4ghz = 3.5 mJy. 

The two-point angular correlation function, w{9), of 
NVSS sources has been measured by Blake & Wall (2002a, 
2002b) and Overzier et al. (2003) for different flux-density 
thresholds between 3 mJy and 500 mJy. The overall shape 
of w{9) is well reproduced by a double power- law. On scales 
below ~ 0.1°, the steeper power-law reflects the distribution 
of the resolved components of single giant radio sources. On 
larger scales the shallower power-law describes the corre- 
lation between distinct radio sources. Since the behaviour 
on small scales is mainly determined by the joint effect of 
the astrophysical properties of the sources and of the reso- 
lution of VLA in the various configurations, while that on 
larger scales is of cosmological origin, we will concentrate 
on the latter only. In fact, the clustering behaviour on scales 
> 0.1° provides insights on both the nature of the radio 
sources, through the way in which they trace the underlying 
dark matter distribution, and on the cosmological frame- 
work which determines the distribution of dark matter at 
each epoch. 

We will consider the two-point angular correlation func- 
tion as measured by Blake & Wall (2002b) for sources with 
SiAGHz > 10 mJy. Such a flux limit ensures the survey to 
be complete and, at the same time, provides good enough 
statistics to measure angular clustering up to very large 
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Figure 1. Adopted redshift distribution per unit redshift inter- 
val, -Sr^z), for the radio source population with 5i.4ghz > 10 mjy. 
The dashed line represents the contribution of AGN-powered ra- 
dio sources according to the PLE model of DP90, while the dotted 
line shows the contribution from star-forming galaxies obtained 
by Magliocchetti et al. (2002). 



scales with small uncertainties. Moreover, this limit also rep- 
resents the deepest flux-density for which systematic sur- 
face density gradients are approximately negligible (Blake 
& Wall 2002a). The NVSS source surface density at this 
threshold is 16.9 deg~^. 

The data corresponding to separations in the range 
0-1° i$ ^ ^ 0.3° most likely suffer from a deficit of pairs 
probably due to an imperfect cleaning of bright side lobes 
(sec Blake, Mauch & Sadler, 2004). Therefore in the follow- 
ing we will only consider scales 9 > 0.3°. In this range of 
scales, the measured angular correlation function is found 
to be described as a power law, w{6) = 1.49 • 10~^ x d~^ °^, 
with 6 in degrees (see Blake, Mauch & Sadler 2004). 



3 REDSHIFT DISTRIBUTION OF 
MILLIJANSKY RADIO SOURCES 

In order to provide theoretical predictions for the angular 
two-point correlation function of a given class of objects it 
is necessary to know their redshift distribution, Af{z), i.e. 
the number of objects per unit comoving volume as a func- 
tion of redshift. Unfortunately, the redshift distribution of 
mJy radio sources is not yet accurately known as the ma- 
jority of radio sources powered by Active Galactic Nuclei 
(AGN) - which dominate the mJy population - axe located 
at cosmological distances (< z >~ 1) and are in general 
hosted by galaxies which are optically extremely faint. 

A large set of models for the epoch-dependent Radio 
Luminosity Function (hereafter RLF) arc available in liter- 
ature (see e.g. Dunlop & Peacock 1990, Rowan-Robinson et 
al. 1993, Toffolatti et al. 1998, Jackson & Wall 1999, Willot 
et al. 2001), but they all suffer from the fact that they are 
mainly based on, and constrained by data sets which include 
only bright sources {Siaghz ^ 100 mJy) so that any extrap- 
olation of their predictions to lower fluxes is quite uncertain. 
Amongst all the available RLFs, those provided by Dunlop 
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& Peacock (1990, hereafter DP90) are the most commonly 
used to infer the redshift distribution of radio sources at the 
mJy level. DP90 derived their set of RLFs on the basis of 
spectroscopically complete samples from several radio sur- 
veys at different frequencies. By using a 'maximum entropy' 
analysis they determined polynomial approximations to the 
luminosity function and its evolution with redshift which 
were all consistent with the data available at that time. In 
addition, they also proposed two models of a more physical 
nature, assuming cither Pure Luminosity Evolution (here- 
after PLE) or Luminosity /Density Evolution. 

Recently, with the aim of determining the photometric 
and spectroscopic properties for at least the population of 
local (i.e. z < 0.2) radio sources at the mJy level, a rmmber 
of studies have concentrated on samples obtained combining 
radio catalogues like FIRST and NVSS with optical ones like 
SDSS and 2dFGRS (Sadler et al. 2002, Ivezic et al. 2002, 
Magliocchetti et al. 2002, 2004). The radio/optical samples 
obtained in this way provide a crucial constraint at z ^ 
for any theoretical model aiming at describing the epoch- 
dependent RLF. 

For instance, Magliocchetti et al. (2002, hereafter M02) 
have shown that the RLF at 1.4 GHz derived from all the 
objects in their spectroscopic sample having 5 > 1 mJy 
and bj < 19.45, is well reproduced at relatively high radio 
luminosities (Pi.4GHz > 10^°-^ WKz"! sr"^) by the DP90's 
PLE model for steep-spectrum FRI-FRII sources (Fanaroff 
& Riley 1974). At lower radio luminosities, where the ra- 
dio population is dominated by star-forming galEixies, the 
measured RLF is better described by the one proposed by 
Saunders et al. (1990) for IRAS galaxies (see also Rowan- 
Robinson et al. 1993), although with a small adjustment of 
the parameters. 

Following to M02, we adopt the DP90 PLE model to de- 
rive the redshift distribution of AGN-fuelled radio sources 
with S'i.4GHz > lOmJy (see Fig. 1, dashed line). We will 
not take into account the other DP90 models, since we 
have found them to be inconsistent with the local RLF. In 
fact, while the model assuming density/luminosity evolution 
over-predicts the number of steep-spectrum radio galaxies 
below Pl.4GHz < 10^^ WHz-isr-\ models using polynomial 
approximations for the RLF give an unrealistic overestimate 
of the number of local sources at every luminosity. 

We will also use the fit provided by M02 to the local 
RLF in order to estimate the redshift distribution of star- 
forming radio galaxies with S'i.4ghz > lOmJy (see Fig. 1, 
dotted line). Note that star- forming galaxies comprise less 
than 30% of the « 5 0.1 population of radio sources with 
iSi.4GHz > lOmJy, and only ~ 0.5% of the total counts at 
this flux limit. As we will show in the next Section, this re- 
sult implies that the contribution of star-forming galaxies to 
the large scale angular clustering of NVSS sources is totally 
negligible. 



4 MODEL FOR THE ANGULAR 
CORRELATION FUNCTION 

The two-point angular correlation function, w{6), of a popu- 
lation of extragalactic sources is related to their spatial cor- 
relation function, ^{r,z), and to their redshift distribution 
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via the Limber's (1953) equation: 

w(6>) = / dzAf^iz) / d{5z)^[r{6z,e),z] 

Jz Jz'(z) 



dzAfiz) 



(1) 



In this expression, r{Sz,6) is the comoving spatial distance 
between two objects located at rcdshifts z and z + Sz and 
separated by an angle 9 on the sky. For a flat universe and 
in the small angle approximation (which is still reasonably 
accurate for the scales of interest here, i.e. 0.3° ^ 6 ^ 10°), 



Hiz) 



5z^ +dl{z), 



(2) 



where H{z) = HoE{z) is the time-dependent Hubble pa- 
rameter and de(«) is the comoving linear distance on the 
sky surface corresponding, at a the redshift a, to an angular 

separation 6. Integrations in Eq. (1) are performed in the 

ranges 2 = [Zmin, Z^a.^] and Z'{z) = [Zmin - ^max - z], 

where Zmin = and Zmax ~ 6 are, respectively, the mini- 
mum and the maximum redshift at which radio sources are 
observed. 

On sufHcicutly large scales (e.g. r ^ 3Mpc), where the 
clustering signal is produced by galaxies residing in distinct 
dark matter halos and under the assumption of a one-to- 
one correspondence between sources and their host halos, 
the spatial two-point correlation function can be written as 
the product of the correlation function of dark matter, ^dm, 
times the square of the bias parameter, b (Matarrese et al. 
1997, Moscardini et al. 1998): 



^{r,z) = b%M^s,z)^T,Mir,z). 



(3) 



Here, M^g represents the effective mass of the dark matter 
haloes in which the sources reside and b is derived in the 
extended Press & Schechter (1974) formalism according to 
the prescriptions of Sheth & Tormen (1999). 

The function ^dm is determined by the power spectrum 
of primordial density perturbations as well as by the under- 
lying cosmology. For the power spectrum of the primordial 
fluctuations wc adopt the fltting relations by Eisenstein & 
Hu (1998) which account for the effects of baryons on the 
matter transfer function. The initial power spectrum is as- 
sumed to be scale invariant with a slope n = 1. As already 
mentioned, we adopt a 'concordance cosmology', in agree- 
ment with the flrst-year WMAP results (Spcrgcl et al. 2003). 

In the range of scales of interest here the cluster- 
ing evolves in the linear regime, so that ^DM{r,z) = 
D^iz)i,DM{r,0), D{z) being the linear density growth-rate 
(Carroll, Press & Turner 1992). 

Since the population of radio sources with Si.4ghz > 10 
mjy is composed by two difl'erent types of objects, i.e. AGNs 
and star-forming galaxies, which display difiierent clustering 
properties (see e.g. Saunders, Rowan-Robinson & Lawrence 
1992, Madgwick et al. 2003), the angular correlation func- 
tion of the whole NVSS sample is given by [cf. e.g. Wilman 
et al. 2003, Eq. (9)]: 



w{6) = /agnWagn(^) + /sfWsf(^) + 2/agn/sfWct 



.(^) (4) 



where /agn and /sf are the fractions of AGNs and star- 
forming galaxies in the whole sample, respectively, i.e. 



J^dzJ\fAGN/SF{z) 

Jagn/sf — 



dzJ\f{z) 



(5) 



Af{z) being the global redshift distribution. 

In Eq. (4), wagn and wsf are the angular correlation 
functions of the two classes of radio sources, while uicross 
accounts for the cross-correlation between the two popula- 
tions: 



-L 



dzAfAGN{z)AfsF{z) 



L 



d{5z)(, cross [r{5z,e),z\ 



We model 
Eq. (31)]: 



Z'(z) 

j dzMAG¥i{z) j dzJ\fsF{z 

J Z J Z 

as [cf. e.g. Magliocchetti et al. 1999 



(6) 



^{r,z) = \/CAGN(r,2;)^sF(r,2;) 



biMti"*, z)b{M^i, z)^DMir, z) 



(7) 



where ^agn and ^sf are the spatial correlation functions of 
AGNs and star-forming galaxies respectively, while M^g^^ 
and M^g denote the effective masses of the corresponding 
dark matter halos [cf. Eq. (3)]. Note that this deflnition 
likely overestimates the cross-correlation term, which may 
be close to zero, since AGN-powered radio sources are nor- 
mally found in clusters while star-forming galaxies preferen- 
tially reside in the field. 



5 RESULTS 

As a first step, we obtain estimates of the angular two-point 
correlation function, w{9), assuming a redshift-independent 
effective mass for the host halos of AGN-powered radio 
sources, , consistent with the clustering properties of 

optically selected quasars (Porciani, Magliocchetti & Nor- 
berg 2004; Groom et al. 2005). The value of M^^'^ is taken 
as a free parameter. The observationally determined spatial 
correlation length of star- forming galaxies, r(,{z = 0) ~ 3-4 
Mpc/h (Saunders, Rowan- Robinson & Lawrence 1992; see 
also Wilman et al. 2003) is consistent with an effective halo 
mass not exceeding M^g = 10^^ Mq /h. We adopt this value 
in our analysis. Any dependence of M^g on redshift can be 
ignored because of the small range in redshift covered by 
star-forming galaxies with Si.4GHz > 10 mJy. 

In the top panel of Fig. 2 we show the contributions 
to w{6) of NVSS sources with S'i.4ghz > lOmJy arising 
from the clustering of AGNs with M^^^ = 10^^, lO^"*, and 
10^ ' Mq/Ii (short-dashed curves) and star-forming galaxies 
(long-dashed curve), as well as from the cross-correlation 
between these two populations (dotted curves). As noted 
above, the cross correlation term, computed from Eq. (7) 
for the three values of AI^^ , is likely overestimated. The 
above contributions are then compared with the observa- 
tional determination of Blake & Wall (2002b). In the lower 
panel of the same Figure the redshift-evolution of the bias 
parameter is represented. 

Figure 2 shows that the contribution of star-forming 
galaxies to the observed angular correlation function is neg- 
ligible. This is because their contribution to the S'i.4GHz > 10 
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Figure 2. Top panel: contributions to the two-point angular 
correlation function of NVSS sources with Si.4ghz > 10 mjy 
from AGNs (sliort-dashed curves), star-forming galaxies (long- 
dashcd curve), and their cross-correlation (dotted curves). The 
effective halo mass of AGN-powered radio galaxies is assumed 
to be constant with redshift. The three short-dashed and dotted 
curves correspond, from bottom to top, to M^^'^ = 10^^, 10^'' 
and 10^^ Mq /h, respectively. The elfective mass of the halo in 
which star-forming galaxies reside is = 10^^ Mq/Yi (long- 

dashed curve). Points with error bars represent the angular two- 
point correlation function of NVSS sources as measured by Blake 
& Wall (2002b). Lower panel: Evolution of the bias parameter 
for AGNs (short-dashed curves) and star-forming galaxies (long- 
dashed curves). 



mJy counts is small, implying /sf ~ 5 x 10~^ [see Eq. (4)], 
compared with /agn ~ 1- We note that this conclusion is in 
disagreement with that of Cress & Kamionkowsky (1998), 
which relies on the adoption of a local RLF of star-forming 
galaxies noticeably exceeding the direct estimates of M02. 

Clearly, models in which the mass of the halo hosting 
AGNs is constant in redshift badly fail at reproducing the 
overall shape of the observed angular correlation function. 
This is because contributions to wagn(9) on a given angular 
scale come from both local and high-redshift sources. For a 
redshift-indopcndcnt halo mass, on scales 6^2° the positive 
contribution of local sources is overcome by the negative con- 
tribution of distant AGNs, which dominate the redshift dis- 
tribution and whose bias parameter rapidly increases with 
redshift [see Eq. (3) and Fig. 2]. In fact, for redshifts a ^ 0.5- 
1 such angular scales correspond to distances where fAGN(r) 
is negative. 

This situation is illustrated by Fig. 3, showing that 
the contribution to wagn(^) from local AGNs {z < 0.2, 
dot-dashed hue) wfth M^'^ = 10^^ M© /h comes close to 
accounting for the observed angular correlation function 
(the dot-dashed line also includes the contribution of star- 
forming galaxies with M^g = lO^^Mo/h, which however 
turns out to be almost negligible). But if the same value of 
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Figure 3. Contributions to the two-point angular correlation 
function of NVSS sources with Si.4ghz > 10 mJy from radio 
sources above (dashed and dotted lines) and below (dot-dashed 
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Star-forming galaxies contribute (yet marginally) only at z < 0.2. 
The contribution of higher redshift sources (AGNs only) is posi- 
tive for (? < 2° , and becomes negative on larger scales (the dotted 
line shows its absolute values). 



also applies to high-redshift AGNs, wagn(6') turns 
out to be too high on small scales and negative on large 
scales. Therefore, the only way to account for the behaviour 
of the data on large scales is to weight down the high-« con- 
tribution by decreasing the bias factor at such redshifts. We 
note that a similar result - which drastically rules out mod- 
els featuring a bias function strongly increasing with z - was 
obtained by Magliocchetti et al. (1999) via their counts in 
cells analysis of the FIRST data. 

Having shown that the contribution of star-forming 
galaxies to the clustering signal is negligible, hereafter we 
will only consider AGNs and will drop the 'AGN' index from 
the effective mass. 

Based on the fact that AGN-powered radio galaxies are 
mainly found in very dense environments such as groups 
or clusters of galaxies, we can then assume that the epoch- 
dependent effective mass controlling their clustering proper- 
ties, Meff is proportional to the characteristic mass of viri- 
alized systems, Mi,{z), which decreases with increasing red- 
shift. M^,{z) is defined by the condition (Mo & White 1996): 



a[M4z),z] = 5ciz) 



(8) 



where a{M, z) is the rms fluctuation of the density field at 
redshift z on a. mass scale M, while 5c{z) is the threshold 
density for collapse of a spherical perturbation at redshift z. 

then represents the typical mass scale at which the mat- 
ter density fluctuations collapse to form bound structures. 
As required by the data, this choice for McS (whose redshift 
dcpcudcucc is shown in the bottom panel of Fig. 4) leads to 
an almost redshift-independent bias factor (central panel of 
Fig. 4) except for the highest value of Meff(O) for which we 
have a moderate decrease of b with increasing z. 

This model successfully reproduces the observed w{d) 
up to scales of at least ~ 4° (top panel of Fig. 4) . The best 



fit is obtained for Meff(z = 0) = 10^ 



Mq/Yi. On the 



largest angular scales {9 ^ 5°), the model correlation func- 
tion still falls short of the data. However, small systematic 
variations in the source surface density due to calibration 
problems at low fiux densities may spuriously enhance the 
estimate of w{9) (Blake & Wall 2002a). An offset as small 
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Figure 4. Two-point angular correlation function (top panel), 
bias parameter (central panel), and effective mass of the host 
dark matter halo (bottom panel) for radio sources with S'i.4ghz ^ 
lOmJy for a model in which M^ff ex Mi,{z) (see text for details). 
The three dashed curves correspond, from bottom to top, to dif- 
ferent local values of the effective mass: Meff(^ = 0) = 10'3, 10'4 
and 10^^ Mq /h. The solid line in the top panel represents the 
prediction obtained for M^siz = 0) = 10^*-^ Mq/Tci by adding 
a constant offset e = 0.0001 to w{6) in order to account for the 
effect of possible spurious density gradients in the survey. The 
dotted curves represent the prediction for flo,m = 0.105 and 
Meftiz = 0) = IQi^-lSMo/h. Points with error bars show the 
two-point angular correlation function of NVSS sources as mea- 
sured by Blake & Wall (2002b). 



as £ = 10 would be enough to remove the discrepancy 
between the model and the data (solid line in the top panel 
of Fig. 4). 

On the other hand, if the measured w{9) is unaffected by 
spurious density gradients in the survey, it may be possible 
to reconcile the model with the observed signal by exploiting 
the dependence of the model predictions on cosmological 
parameters, and notably on the matter density parameter, 
on the baryonic fraction and on the Hubble constant (see 
Matsubara 2004, Fig. 1). Decreasing f2o,m shifts the cut- 
off of the dark matter spatial correlation function ^dm to 
larger scales, while a non-null value of Qb determines a peak 



in the correlation function just below the cut-off scale. The 
amplitude of this peak is an increasing function of Qb- 

A change in the value of h has an effect which is qual- 
itatively similar to that induced by changing f2o,m, but of 
a significantly smaller amplitude; for this reason we have 
decided to keep this parameter fixed to the reference value 
h = 0.72. We have also ignored the dependence of the clus- 
tering pattern on the rms mass fluctuation, as' in principle 
this quantity contributes to the overall normalization of the 
spatial clustering, but we have checked that in our case the 
effect of changing it is almost negligible. This is because 
the decrease of b with increasing as is compensated by the 
corresponding increase of the effective mass (assumed to be 
proportional to l\/h{z)). Therefore we also keep as = 0.9. 

The parameters in our analysis are then: Meff(O), Qo.m 
and fib- In Fig. 5 we show the predicted two-point angular 
correlation function for different choices of these parameters. 
In each panel, the dashed curves correspond (from bottom 
to top) to Qb = 0.03, 0.045 and 0.06. The values adopted 
for both MefF(O), and Oo,m are given in each panel: Mo£f(0) 
increases from the top to the bottom panels, while Qo,m 
increases from the left-hand to the right-hand panels. 

We note that variations of these parameters within the 
allowed ranges mainly translate into changes of the normal- 
ization of w{6), although varying i^b also affects the w{9) 
slope on scales 6^2°, while changes in f2o,m modify its slope 
on larger scales. In general, the amplitude of w{6) increases 
with increasing Afoff (0) or fit,, and decreases for decreasing 
f2o,m- From Fig. 5 we can deduce that, given the observed 
w{6), smaller values of Mgff (O) are favoured for lower values 
of Ho,™. 

Fixing the baryonic content to the reference value Q.b = 
0.045, we have then investigated the Ocm-Moff (O) interde- 
pendence by constructing contours on the fio,m-Afefi(0) 
plane. The results are shown in Fig. 6 where the curves rep- 
resent, from the innermost to the outermost, contours cor- 
responding to the 68.3, 95.4 and 99.7 per cent confidence 
intervals, respectively. The filled square corresponds to the 
best-fit set of parameters: no,^ = 0.105, Mee{z = 0) = 
-[^q13.15 jy^^^j^ The resulting angular two-point correlation 
function, and the redshift-evolution of both the bias param- 
eter and the effective mass are represented by the dotted line 
in Fig. 4. We note that the effect of a change in the baryonic 
content is just that of smearing the relation between fio.m 
and Mcff (O) (in particular for f2o,m ~ 0.2) since there is a 
whole set of fio,m-Afeff(0) pairs which can provide the same 
best-fit to the data. 

The fit obtained for a zero offset e is really good, but 
is obtained for a local matter density parameter Oo.m sub- 
stantially smaller than what indicated by other data sets 
(see e.g. Spergel et al. 2003), although the difference is sig- 
nificant to less than 2>a. Taken at face value, our results 
indicate that, within the standard cosmological framework, 
in the local universe radio sources arc as strongly clustered 
as rich clusters of galaxies, in excellent agreement with the 
findings of Peacock & Nicholson (1991) and Magliocchetti 
et al. (2004). 

As a last comment, we note that the degeneracy be- 
tween the host halo mass and the matter density parameter 
indicates that large area radio surveys without redshift infor- 
mation can hardly provide significant constraints on cosmo- 
logical parameters, unless independent information on Mefi 



© 0000 RAS, MNRAS 000, 000-000 



The large scale clustering of radio sources 7 



0.01 

g 0.001 

0.0001 

10-5 
0.01 

g 0.001 

0.0001 

10-s 
0.01 

g 0.001 

0.0001 



10 



-5 



I I I 1 1 1 1 I I I I I I I 



M,„(0)=10'3: 
no,„,=o.n 







\- 


















— h 



Me„(0)=101' 























— 










— 



iii. -\ 



I Mill 



J I I I I I I 



I I I I I I I I I I I I I I : 

M,„(0)=10'3: 



U,m 



'hi 



Me„(0)=10" 



\1 



M,„(0)=10' 



I Mill 



J I I I I I I I 



I I I I I I I I I I I I I I I = 

M,„(0)=10'3: 



n„^=o.3" 



in 



-ml — I t i I N I 



M,„(0)=10i^= 



-m\ — I i lj iiii| 



M,„(0)=10'=^ 



o,m 



I Mill 



..... 

I I I I I- 



0.5 1 5 10 

0(deg) 



0.5 1 5 10 

0{deg) 



0.5 1 5 10 

0{deg) 



Figure 5. Dependence on cosmological parameters of the two-point angular correlation function of radio sources with 5i.4ghz > 10 
mjy. The effective mass of the dark matter haloes in which the sources reside is assumed to be proportional to M*(z). In each panel, 
the adopted values of Mo£f(0) (in solar units) and of the matter density parameter, f2o,m are indicated: the former increases from top 
to bottom panels while the latter increases from left-hand to right-hand panels. In each panel the three dashed curves correspond to 
different choices for the baryonic mass density parameter: Qi, = 0.03, 0.045 and 0.06 from the lower to the upper curves, respectively. In 
the bottom right-hand panel the dotted curves (almost indistinguishable from each other) represent the two-point correlation function 
obtained by adding to the dashed curves a constant offset e = 10~* (see text for details). 



is available. On the other hand, for a given set of cosmolog- 
ical parameters, such surveys yield important information 
on the masses of the hosting halos. 



6 DISCUSSION AND CONCLUSIONS 

The observed angular two-point correlation function of mJy 
radio sources exhibits the puzzling feature of a power-law 
behaviour up to very large (~ 10°) angular scales. Stan- 
dard models for clustering, which successfully account for 
the clustering properties of optically selected quasars, turn 
out to be unable to explain the long positive tail of the 
w{6), even when 'extreme' values for the parameters are in- 
voked. This is because - according to the standard scenario 



for biased galaxy formation in which extragalactic sources 
are more strongly clustered at higher redshifts - the clus- 
tering signal of radio sources at 2; > 1, which is negative on 
large scales, overwhelms that of more local sources, yielding 
a sharp cut-off in the angular two-point correlation function 
on scales of ~ 1-2 degrees. 

The only way out we could find is to invoke the cluster- 
ing strength of radio sources to be weaker in the past. The 
data can then be accounted for if we assume that the char- 
acteristic mass of the halos in which these objects reside, 
Meg{z), is proportional to M*(«), the typical mass-scale at 
which the matter density fluctuations collapse to form bound 
structures (see Mo & White 1996). 

A good fit of the observed w{6) up to scales of at least 
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4° degrees is obtained for M^h{z = 0) = io^*-^'^o °'^ Mo/h. 
The data on larger scales can be accurately reproduced if 
the measured values of w{d) are slightly enhanced by small 
systematic variations in the source surface density due to 
calibration problems at low flux densities (Blake & Wall 
2002a). In the absence of such a systematic offset, the data 
might indicate a lower value of the mean cosmic matter den- 
sity than indicated by other data sets. The best fit is indeed 
obtained for Qo.m = 0.105 and Mes{z = 0) = lO^^-^^ Me/h, 
with rather large uncertainties on both parameters, as shown 
by Fig. 6. In particular, the best fit value of Qo.m is less than 
3(7 away from the 'concordance' value. This shows that cur- 
rent large scale radio surveys without redshift measurements 
cannot provide strong constraints on cosmological parame- 
ters, but are informative on the evolution history of the dark 
matter halos hosting radio sources. 

Taken at face value, the above results point to different 
evolutionaxy properties for different types of AGNs as the 
decreasing trend for the effective mass ruling the clustering 
of radio sources found in this work strongly differs from the 
behaviour of that associated to optically selected quasars. 
For the latter sources, the data are in fact consistent with 
a constant Mgg ~ lO'^^ Mq /h, at least up to the highest 
probed redshift, z ~ 2.5 (see e.g. Grazian et al. 2004; Por- 
ciani et al. 2004; Groom et al. 2005) This different behaviour 
is illustrated in Fig. 7, where the values of the bias inferred 
for optically selected quasars have been scaled to the values 
of as and flo,m used here. 

It must be noted, however, that the present data provide 
very weak constraints on the clustering properties of radio 
sources at 2; > 1. For example, we have directly checked 
that the predicted w{0) does not change significantly over 
the range of angular scales considered here if we assume 
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Figure 7. Bias factor as a function of redshift. The data points 
are from observations of optically selected quasars (open squares 
from Porciani et al. 2004, filled dots from Croom et al. 2005, 

open triangle from Grazian et al. 2004; the values of b derived 
by these authors have been scaled to the values of as and Qo^rn 
used here), while the dashed line represents the functional form 
derived in this work for the mjy population of radio sources. The 
dotted curve shows the bias function obtained for a characteristic 
halo mass constant in time and equal to ~ 10^'^ Mg/h. 



Afoff oc Af*(z) for z < 1.5 and M^a ~ const ~ 10^^ A^©/h at 
higher redshifts. It is thus possible that the difference in the 
clustering evolution between AGN-powered radio sources 
and optically selected quasars is limited to z < 1, consistent 
with observational indications that, at higher redshifts, the 
environment of radio-quiet and radio-loud QSOs is almost 
the same (Russel, Ellison & Benn 2005). 

Interestingly, our analysis indicates that, at « ~ 1, the 
effective halo mass associated to radio sources is consis- 
tent with being essentially equal to that of optical quasars 
(Croom et al. 2005). This suggests that, at least in this red- 
shift range, the bias parameter for radio-loud and radio-quiet 
AGNs is similar. 

On the other hand, our findings suggest that, at least for 
z < 1 and at variance with what found for optical quasars, 
the clustering of radio sources refiects that of the largest 
halos which collapse at any given cosmic epoch. This con- 
clusion is in keeping with results of previous studies showing 
that, locally, radio sources arc preferentially associated with 
groups and clusters of galaxies (e.g. Peacock & Nicholson 
1991; Magliocchetti et al. 2004), and, at higher redshift, are 
often associated with very massive galaxies and very massive 
galaxy environments (e.g. Carilli et al. 1997; Best, Longair 
& Rottgering 1998; Veriemans et al. 2002; Best et al. 2003; 
Groft et al. 2005; Overzier et al. 2006). 

The major intriguing point that remains unsolved is 
the link with the population of optical QSOs. The cluster- 
ing properties of these latter objects seem to reflect that of 
"normal" elliptical galaxies. So why this difference? Glearly, 
more observations are crucial to shed light on this issue. So 
far the main limitation to our understanding of the environ- 
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mental properties of radio sources has been due to selection 
effects that allow identifications of mostly radio galaxies in 
the local universe and exclusively quasars (regardless of their 
radio activity) at higher redshifts. Therefore, it would be of 
uttermost importance to consider a redshift range in which 
the clustering properties of both radio galax;ies and radio- 
active quasars can be measured with good precision. We 
plan to tackle this issue in a forthcoming paper. 
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